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Share price returns on different time scales can be well modeled by a superstatistical dynamics. 

Here we provide an investigation which type of superstatistics is most suitable to properly describe 
share price dynamics on various time scales. It is shown that while x^-superstatistics works well on 
a time scale of days, on a much smaller time scale of minutes the price changes are better described 
by lognormal superstatistics. The system dynamics thus exhibits a transition from lognormal to x^ 
superstatistics as a function of time scale. We discuss a more general model interpolating between 
both statistics which fits the observed data very well. We also present results on correlation functions 
of the extracted superstatistical volatility parameter, which exhibits exponential decay for returns 
on large time scales, whereas for returns on small time scales there are long-range correlations and 
power-law decay. 


I. INTRODUCTION 

Many well established concepts in mathematical fi¬ 
nance (such as the Black-Scholes model) are based on 
the assumption that an index or a stock price follows a 
geometric Brownian motion, and as consequence the log 
returns of these processes are Gaussian distributed. But 
nowadays it is well known that the log returns of realis¬ 
tic stock prices are typically non-Gaussian with fat tails 
Q-dl. Such behaviour can be well captured by super- 
statistical models i-El. The basic idea of this method 
borrowed from nonequilibrium statistical mechanics is to 
regard the time series as a superposition of local Gaussian 
processes weighted with a process of a slowly changing 
variance parameter, often called /3. This approach has 
been applied to many areas of complex systems research, 
including turbulence, high energy scattering processes, 
heterogenous nonequilibrium systems, and econophysics 
(see e.g. [ll| for a short review). In finance early appli¬ 
cations of the superstatistics concept were worked out by 
Duarte Queiros et al. (6,[3 and Ausloos et al. Q. Van 
der Straeten and Beck Q analysed daily closing prices of 
the Dow Jones Industrial Average index (DJI) and the 
SP 500 index. They verified that both log-normal super- 
statistcs and x^ superstatistics result in good approxima¬ 
tions. Biro and Rosenfeld 0 also studied the data sets 
of the Dow Jones index and verified that the distribu¬ 
tion of log returns is well fitted by a Tsallis distribution. 
Katz and Li Tian [l[ showed that the probability distri¬ 
butions of daily leverage returns of 520 North American 
industrial companies during the 2006-2012 financial cri¬ 
sis comply with the g-Gaussian distribution which can be 
generated by x^ superstatistics. They also verified in 0 
that the Tsallis entropic parameter q obtained by direct 
fitting to g-Gaussians coincides with the q obtained from 
the shape parameters of the x^ distribution fitted to the 
histogram of the volatility of the returns. Gerig, Vicente 
and Puentes Q consider a similar model that indicates 
that the volatility of intra day returns is well described 
by the x^ distribution, see also Q for related work in this 
direction. 

In this paper, we will carefully analyse for various data 


sets of historical share prices which type of superstatis¬ 
tics is best suited to model the dynamics. While Tsallis 
statistics (= g-statistics) is known to be equivalent to x^ 
superstatistics [H, [l^, there are other types of super¬ 
statistics, such as lognormal superstatistics and inverse 
X^ superstatistics [^, which are known to be different 
from g-statistics (though all these different statistics gen¬ 
erate similar distributions if the variance of the fluctua¬ 
tions in P is small El)- We show that in our analy¬ 
sis x^-superstatistics appears best suitable to describe 
the daily price changes, whereas on much smaller time 
scales of minutes lognormal superstatistics seems prefer¬ 
able. We analyse the relevant time scale of the changes 
in the superstatistical parameter /3 and present results 
for the decay of correlations in /3. For small return time 
scales, correlation functions exhibit power law decay and 
there are long memory effects. In the final section, we 
develop a synthetic stochastic model that fits the data 
well. This is kind of a hybrid model interpolating be¬ 
tween lognormal and x^-superstatistics. 

This paper is organized as follows. In section II we 
look at share price returns on large (daily) time scales. In 
section III we do a similar analysis on small (minute) time 
scales. In section IV we investigate correlations of the 
superstatistical volatility parameter on both time scales. 
In section V the hybrid model is introduced. Our final 
concluding remarks are given in section VI. 

II. SUPERSTATISTICS OF LOG-RETURNS OF 

SHARE PRICES ON A LARGE TIME SCALE 

Non-equilibrium system dynamics can often be re¬ 
garded as as superposition of a local equilibrium dynam¬ 
ics and a slowly fluctuating process of some variance vari¬ 
able P [l^. These types of ‘superstatistical’ nonequilib¬ 
rium models are also useful for financial time series [a, l7| . 
In this article, the empirical data we use as an example 
is the historical stock prices of Alcoa Inc(AA), which is 
an American company that engages in the production 
and management of primary aluminium, fabricated alu¬ 
minium and alumina. We have looked at shares of many 
other companies as well (see Tab. 1 in section IV), with 
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similar results. Our data set covers the period January 
1998 to May 2013. We study the log return Ri denoted 
by 

(^) 

where i = 0,1, 2,..., TV; Si and 5^+1 are two successive 
daily closing prices. We consider the normalised log re¬ 
turns 




( 2 ) 


which have been rescaled to have variance 1. The symbol 
(• • •) denotes the long-time average. 

From the simplest superstatistics model point of view, 
the entire time series of stock prices can be divided into 
n smaller time slices T. We call T optimal window size. 
Within each T, the financial volatility /3 is temporarily 
constant and the log return of the stock price is Gaussian 
distributed. /3 has some probability distribution / (/3) to 
take a particular value in a given slice. The conditional 
probability p{u\j3) is 


p{u\^) = 



( 3 ) 


and the marginal probability distribution of u for long 
time observation is the average over local Gaussians 
weighted with the probability density /(/3) 


p{u) = J p{u\l3)f{l3)dl3. 


( 4 ) 


The integration over /3 yields non-Gaussian behaviour 
with fat tails. 

We now describe our technique to obtain the optimal 
window size T for a given time series. Firstly we split 
the time series into 


n = 


[—J 


( 5 ) 


equal intervals, where [ J denotes the floor function and 
At is the dimensionless window size, i.e. the number of 
data points in a given window. N is the total number 
of data points of the entire time series. Generally the 
kurtosis of a random variable u is defined as 


K = 


(u2)2 


( 6 ) 


and it is equal to 3 for a Gaussian distribution of arbitrary 
variance. For a given window size At, the kurtosis in the 
jth window is given by 


KAtU) 


1 4 

AiA,i={j-l)At+l '^i 


( 1 

At — 1 ) 



( 7 ) 
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FIG. 1: Determination of the optimal window size for the 
Alcoa share price data. The intersection with the line kurtosis 
R, — 3 yields T = 18 ± 0.5. The various values of R for a 
given At (indicated by different colors in the online version) 
are obtained for different translational shifts of the sliding 
windows. The scattering of the data can be used to estimate 
the standard deviation as 5r ~ 0.03. 


where j = 1, 2,..., n. When we have all the values of kur¬ 
tosis for all windows, we can calculate an average kurtosis 
of the n windows as 

1 " 

i^At = - (8) 

n 


The aim is to achieve an optimum window size such 
that for a given data set the distribution in each window 
is as close as possible to a Gaussian, but with varying 
variance. For this purpose the optimal window size T 
should satisfy the condition 


RAt — 3. 


( 9 ) 


Fig.[T]shows how the average kurtosis changes with the 
window size. We obtain from condition 0 the optimal 
window size 18 ± 0.5 for this example. The result makes 
financial sense. 18 trading days correspond to a time 
scale of about 3-4 weeks. It is a typical time scale where 
market volatility changes, due to events such as changes 
in the confidence in the future economic development, 
anticipated interest changes, and so on. See also for 
related work. 

With the given optimal window size, we can now calcu¬ 
late the local volatility parameter /3 in each time interval 
as 


Pk 


1 


1 

T-1 2^i={k-l)T+l 


(Ui 


Uif 


( 10 ) 


where fc = 1, 2,..., n. Note that the variance of u in each 
window is P~^. One can then plot a histogram of the Pk 
and fit it with some suitable model distribution. 

Here we will consider three distributions to be com¬ 
pared with our experimental distribution of P, which 
were previously advocated in [T^ . The first one is the 
-distribution for which f(P) is given by 


fiiP) 




di/2 

^di/2-lg-di/3/2/3o 


( 11 ) 
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FIG. 2: Best possible fits that can be achieved for the distri¬ 
bution of the volatility /? of Alcoa shares (plotted by dots), 
in a log-linear (top) and double logarithmic plot (bottom). 
Blue: distribution fi{P) with di = 1.51,/3o = 2.19, Green: 

inverse distribution /2(/3) with d 2 = 0.45, (Jo = 2.19, Red: 
lognormal distribution /3(/3) with s = 0.87,/r = 0.45,. 


The second one is the inverse x^"distribution where 





d2/2 

^-d2/2-2g-d2/3o/2/3 


( 12 ) 


The third distribution that will be tested is the log¬ 
normal distribution for which the probability density 
function is given by 


p(u) 



FIG. 3: Gomparison of the histogram of u (plotted by dots) 
with the 3 types of superstatistics, integrated with the same 
parameters as in Fig. [2] Blue: x^ Superstatistics pi(u), 
Green: inverse x^ Superstatistics P 2 {u), Red: lognormal Su¬ 
perstatistics P 3 {u). 


p(u) 



FIG. 4: Amended Superstatistics Blue: x^ Superstatistics 
Pi{u) with di — 1.51, do = 2.19, Green: inverse x^ Super¬ 
statistics P 2 (u) with d 2 = 1.2, do = 2.19, Red: lognormal 
Superstatistics P 3 {u) with s = 1.2, p, — 0.65. 


Still for consistency we also need to check the validity 
of Eq. (HI). We thus also compare the original histogram 
of returns u with the following integrals where the pa¬ 
rameters take the same values as in Figl^ 


h (/3) = 


1 

\/^sj3 


exp 


( - 

[ 2 .^ ) 


(13) 


where 

/r = ln/3o-y- (14) 

The do in Eq- (HD), (HH), (fl^ is the mean value of fp, 
given by 


1 

/3o = (d) = -E^fc’ (15) 

^ k=l 


and di,d 2 ,s are parameters. Lognormal superstatistics 
often occurs for complex systems described by a cascad¬ 
ing dynamics [l^ , whereas nnd inverse superstatis¬ 
tics are more common for additive degrees of freedom 
contributing to a fluctuating temperature or inverse tem¬ 
perature flOj . 

We have fitted our experimental histograms /(d) with 
the above distributions. Given do) we vary di,d 2 and s 
of Eq. (fTTl) . (IT^ . (IT^ in order to obtain the optimum 
fit to our observed /(d)- It can be seen in Eig. [2] that 
lognormal, x^" and inverse x^ superstatistics all yield a 
more or less decent fit, though inverse x^“Superstatistics 
seems less favorable. 


Pi{u) = J Y^exp ^-idM^^/i(d)'ld i = l,2,3 

(16) 

As shown in Eig. |3l for the integrated densities x^ su¬ 
perstatistics seems to fit better to the probability density 
of u compared with lognormal superstatistics and inverse 
X^ superstatistics. 

Thus, if independent variation of the volatility param¬ 
eter in each interval is assumed, then the data clearly 
point to x^ superstatistics, equivalent to Tsallis statistics 
[22j| . On the other hand, independence of dfe may not al¬ 
ways be a good approximation. There can be strong cor¬ 
relations of the volatility parameter dfci and variations 
of the time scales where it is approximately constant. 
In that case more complicated dynamics arise, and one 
could then possibly get a better fit for the integrated 
distributions p{u) if other effective parameters are used. 
For this reason, we also allowed the fitting parameters for 
Pi{u),P2{u),P3{u) to take on other possible values. The 
result of this ‘amended superstatistics’ is shown in Fig.m 

After the adjustment, we find in FiglUthat in fact all 
three superstatistics can describe p{u) quite well. To dis¬ 
tinguish between them, one would need much more data 
so that the tail behaviour would be clearer. In practice, 
more data are available if one considers price changes on 
much smaller time scales than days. This will be done in 
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FIG. 5: Determination of the optimum window size for 
the 1-minute data of Alcoa. The intersection with the line 
kurtosis = 3 yields T k. 11. 

the next section. 


III. SHORT TIME SCALES 

Let us extend our analysis to returns on much smaller 
time scales. A change of statistics as a function of the 
time scale considered is a common phenomenon for many 
complex systems, see e.g. [illii for work in this direc¬ 
tion. Hence it is interesting to also consider return data 
on much smaller time scales (say, minutes), and see what 
is similar and what is different as compared to the anal¬ 
ysis of the previous section. Let Si be the stock price for 
every recorded minute, in our example chosen as that of 
Alcoa Inc(AA). The total number of data points is about 
1.5 million. We look at the returns 

(^) 

where r is an integer in units of minutes. The log returns 
are again normalized to variance 1: 



There is one small technical problem for these types of 
data, as the returns are not given overnight but only dur¬ 
ing normal working hours. This can lead to big overnight 
jumps and affect the analysis. For this reason, if Si+r and 
Si are from two successive trading days, we removed the 
corresponding log . t = 1 means the log return 

is extracted every minute. Again we determined the op¬ 
timal window size, using the same technique as in the 
previous section. We obtain T fv 11 (see Fig. 0). 

Again this time scale of about 11 minutes makes sense. 
It is a typical time scale on which new relevant informa¬ 
tion becomes available to the traders, leading to changes 
in the small-scale volatility. It also coincides with typical 
time scales on which observed correlations in short-term 
returns start to decay [IJ. Our results of fitting the three 
types of superstatistics are shown in Fig. [B]151 

As can be seen in Fig. [51 the lognormal distribution is 
by far best fit of /(/?) if the time scale is 1 minute. 


f(/3) 




FIG. 6: Best fits that can be achieved for the distribu¬ 
tion of the short-scale volatility parameter /3 (time scale 
of returns: 1 minute). Blue: distribution fi{l3) with 

di = 0.13,/3o = 6.33, Green: inverse distribution /2(/3) 
with d 2 = 2.83, do = 6.33, Red: lognormal distribution /3(/3) 
with s = 1.11,/r = 1.23, top: log-linear plot, bottom: double 
logarithmic plot. 


p(u) 



FIG. 7: Gomparison of histogram of u (plotted by dots) with 
the integrated superstatistics distributions, using the same 
parameters as in Fig (6] Blue: Superstatistics pi(u), Green: 

inverse y^ Superstatistics P 2 {u), Red: lognormal Superstatis¬ 
tics P 3 {u). None of the curves is a good fit, indicating the 
presence of strong correlations for the volatility parameter 
Pk- 


Fig. 7 shows a clear difference as compared to the daily 
data in Fig. 3: The integrated formula now does not give 
good fits to p{u). The reason is that the /3fc on a time 
scale of minutes are not anymore statistically indepen¬ 
dent, hence random sampling of Gaussians with different 
variance is not appropriate anymore. 

After the free adjustment in the parameters of 
Pi(u),p 2 (u),p 3 (u), again both y^ and lognormal super¬ 
statistics can provide good fits of p{u). See Fig. |8l 

If one does not allow for parameter amendments, then 
we can conclude that there is a transition from to log¬ 
normal superstatistics when the time scale changes from 
1 day to 1 minute. Also, a more general conclusion seems 
to be that the assumption of a sequence of independent 
volatility parameters Pk is not valid, as we are getting in 
general differences between the optimum fit of /(/3) and 
the corresponding fit of p{u) written as an integral over 
Gaussians with the same corresponding parameters. 
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FIG. 8: Amended Superstatistics Blue: Superstatistics 

Pi(u) with di = 0.36, /3o = 6.33, Red: lognormal Super¬ 
statistics P 3 (m) with s = 2.7, fj, = 3.9, Green: inverse x^“ 
superstatistics (d 2 = 0.2,/io = 1.8). 


1.0 
0.8 

C„(() 0.6 
C„(0) 0.4 
0.2 
0.0 

0 10 20 30 40 50 

t 

FIG. 9: Correlation function of log-returns u on a daily time 
scale for AA shares. The time unit of t is days. 



FIG. 11: Correlation function of volatility P for returns on a 
daily time scale for AA shares. 



FIG. 12: Correlation function of volatility /3 for returns on a 
time scale of minutes. 



IV. CORRELATION FUNCTIONS 

For the development of a suitable dynamical model, it 
is very important to look not only at probability densi¬ 
ties but also on correlation functions and memory effects 
In our case there are two types of correlation 
functions: the one of the original data Ui, 

^ N-t 

( 19 ) 

and those of the volatility parameter /3/e, 

.. n—t 

E (20) 

” ^ k=l 

Figs. [MT^ show Cu{t)/Cu{0) and Cp{t)/Cp{Q), both for 
the daily returns as well as for the 1-minute returns. As 
is illustrated in Fig. 9 and 10, Cu{t) decays almost im¬ 
mediately to zero, both for the daily and minute data. 

More interesting is the correlation function Cp{t). We 
did an analysis of the decay rates of correlation functions 
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of the volatility for many different shares from different 
sectors, the results are summarized in Tab. 1. We ob¬ 
serve that the correlation functions of volatility decay 
in an exponential way for daily returns, ^ e”'*'*, 

whereas for minute return there is a power law decay 
(7/3 (f) ^ t~°‘ with a periodic modulation, see Figs. [TTJT2] 
for the example of A A shares. The strongest correlation 
decay (largest 7) on the daily scale is observed for shares 
from basic materials, whereas the power law decay (ex¬ 
ponent a) on the small time scale is largest for healthcare 
shares and shares from the consumer good sector. Note 
that a strong decay of the volatility correlation function 
in a sense measures a ‘volatility of a volatility’ and is an 
interesting quantity to study. The period of oscillations 
that we observe in figures such as Fig. 12 corresponds 
(roughly) to one trading day and is consistent with peri¬ 
odic oscillations of intraday volatility reported previously 


Company 

Sector 

Y 

T [Days] 

a 

T [Mins] 

Alcoa Inc. (AA) 

basic materials 

0.115 

18 

0.094 

11 

The Coca-Cola Company (KO) 

consumer goods 

0.061 

15 

0.101 

13 

Bank of America Corporation (BAC) 

financial 

0.057 

20 

0.073 

13 

Johnson & Johnson (JNJ) 

healthcare 

0.041 

16 

0.102 

13 

General Electric Company (GE) 

industrial goods 

0.068 

18 

0.062 

14 

Wal-Mart Stores Inc. (WMT) 

services 

0.036 

16 

0.092 

13 

Intel Corporation (INTC) 

technology 

0.048 

17 

0.096 

13 


FIG. 10: Correlation function of log-returns u on a time scale Tab. 1: Decay rates of correlation functions of volatility for 
of minutes. The time unit of t is minutes. shares of different sectors 
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V. SYNTHETIC MODEL 

Based on the results of the previous sections, it is de¬ 
sirable to construct a simple superstatistical dynamical 
model that incorporates the possibility of both lognormal 
and superstatistics on different scales, and allows for 
different decay patterns of correlation functions. 

Here we propose the following model. We start from a 
linear superstatistical Langevin equation 

u = —ju + aL{t) (21) 

where Lit) is Gaussian white noise and the ‘inverse 
temperature’ /3, in accordance with Einstein’s theory of 
Brownian motion, is defined as 


Given a fixed /3, the variance {u^{t)) is given by ^{u^) = 
I3~^ for time t —> oo. However, for the superstatisti¬ 
cal version we allow for fluctuations in the parameter /3. 
Then the above Langevin equation is -by construction- 
superstatistical as we do not keep the parameter P con¬ 
stant but regard it as a random variable that fluctuates 
on a large time scale. Eq. m generates a stochastic 
process and in the end, after P has taken on many differ¬ 
ent values, one may rescale the entire time series u{t) to 
variance 1 using the variance of the complete time series, 
as we did in eq.(2) and (18)). 

Let us now consider n -I- 1 Gaussian random variables 
Xi, i = 0,1,2 ,..., n which are statistically independent 
and have the same variance and mean 0 (except for Xq 
which may have potentially a different variance and dif¬ 
ferent mean). We then write P as 

P = T (1 — k)(Y^ -|- X 2 -l- ... -I- X^, (23) 

where n G [0,1] is a parameter. We now see that if 
At = 1, this system generates lognormal superstatistics, as 
log P = Xq is a Gaussian random variable. On the other 
hand, if /t = 0 this system generates x^-superstatistics 
with n degrees of freedom, as in this case P = 
is x^ distributed. Choosing any value of k G [0,1] one 
can interpolate between lognormal and x^ superstatis¬ 
tics, getting a mixed type of behaviour. 

The Gaussian random variables X^ can again be simu¬ 
lated by ordinary linear Langevin equations of the form 

W =-rY, + EL,(t), i = 0,...,n (24) 

For constant T and E these equations generate the Orn- 
stein Uhlenbeck process, i.e. a Gaussian Markov process 
with exponential decay of correlation functions. More 
complicated dynamics, leading e.g. to power law decay 
of correlation functions, can be constructed if the driv¬ 
ing forces in these linear stochastic differential equations 
are not Gaussian white noise but more complicated cor¬ 
related processes, or critical maps with a near-vanishing 
Liapunov exponent 


f(« 



FIG. 13: Mixed distribution fit to f{P) with k = 0.36 on the 
daily time scale. 


f(/3) 



FIG. 14: Mixed distribution fit to f{P) with k = 0.92 on the 
time scale of minutes. 


Fig. [inland Fig. HH show that indeed the observed dis¬ 
tributions of f{P) for Alcoa shares are best fitted by in¬ 
termediate distributions (a superposition of a lognormal 
and x^ distribution with appropriate weights). The pa¬ 
rameter K increases if one goes from larger to smaller time 
scales of returns. The mixed synthetic model is able to 
reproduce the transition scenario of observed densities 
from x^ superstatistics to lognormal superstatistics in a 
quantitatively correct way, giving good fits on any time 
scale. 

We did this analysis for a variety of time scales r of 
returns, taking again the example of Alcoa shares. In 
Fig. [12] we show how the parameter k depends on the 
time scale of returns. As expected, the parameter k that 
best fits the observations decreases as a function of time 
scale. In fact we observe a logarithmic dependence if the 
time scale is not too big, see the straight line fit in Fig. 15. 

One final remark is at order: One may generalize the 
superstatistics concept to more complicated local pro¬ 
cesses that are not locally Gaussian. Indeed, due to cor¬ 
relations present on small time scales, and/or due to a 
lack of clear time scale separation different distributions 
than Gaussians may locally be present. In this case one 
can still superimpose these local distributions by letting a 
suitable variance parameter fluctuate. It is remarkable, 
however, that for the financial data analysed here this 
generalization to more complicated non-Gaussian local 
processes is not necessary: The simplest superstatistics 
model based on local Gaussians fits our data well, assum¬ 
ing the interpolating model eq.(23) where the probability 
distribution of P changes as a function of scale. 
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FIG. 15: Parameter k describing the relative weight of lognor¬ 
mal and superstatistics in the mixed model as a function 
of the time scale r of return, k decreases if the time scale 
T is increased. For not too big time scales r a logarithmic 
dependence is observed: The straight line corresponds to a fit 
of the first six data points of the form k = 0.907 — 0.044 log r. 

VI. CONCLUSION 

Many investigations of complex systems in the past 
have focused on the application of a particular statis¬ 
tics, for example g-statistics [l^, and then studying the 
effect of varying system parameters, which may change 
the entropic index q. Here we have shown that for finan¬ 


cial time series it is sometimes useful to consider broader 
classes of statistics and even proceed from one class of 
superstatistics to another when the scale or other system 
parameters under consideration are changed. The exam¬ 
ple we considered in detail in this paper were share price 
returns of various companies. We provided evidence that 
there is a transition scenario from lognormal superstatis¬ 
tics to superstatistics, with lognormal superstatistics 
giving a better fit to the data on small time scales and 
superstatistics (= g-statistics) on larger time scales. 
We constructed a hybrid superstatistical model that al¬ 
lows to implement both types of superstatistics, with a 
weighting parameter k that describes how far away we 
are from one of the two cases. Correlation functions 
of the extracted superstatistical volatility parameter Pk 
were shown to exhibit different qualitative behavior as 
a function of the time scale of returns, with exponential 
decay on large time scales and power law decay on small 
time scales, modulated by intraday periodicity. The de¬ 
cay parameters of the exponential or power law decay 
were extracted from the data and were shown to depend 
slightly on the sector of shares considered. The general 
transition scenario from lognormal to x^ superstatistics 
as a function of the time scale of returns, however, is a 
general phenomenon and occurs for all sectors in a similar 
way. 
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